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Abstract 



In this work, supersymmetric quantum mechanics wiU be used to obtain complex solutions to 
Painleve IV equation with real parameters. We will also focus on the properties of the associated 
Hamiltonians, i.e. the algebraic structure, the eigenfunctions and the energy spectra. 

^ l' Keywords: quantum mechanics, non-linear differential equations, Painleve equations, complex 

I ' > , potentials with real spectra 

^ ! 1 Introduction 

In recent years a growing interest in studying non-linear models has appeared, since they can be used to 
describe many physical situations. This leads typically to the analysis of non-linear differential equations, 
^ ' which are harder to solve than the linear ones. However, there are solution methods to solve them, 

which mainly consist in reducing non-linear equations to linear ones. This approach has been employed 
for obtaining solutions of the Painleve equations, which we are particularly interested in. Although 



■ discovered from strictly mathematical considerations, the Painleve equations are widely used nowadays 



to describe physical phenomena [T|. 

When is based on supersymmetric quantum mechanics (SUSY QM), this method transforms the non- 
linear equations into Schrodinger type ones. In addition, if we end up with a previously solved Schrodinger 
equation this problem will be further simplified. In particular, this works for the Painleve IV equation 
(Piv) [HIS] J which is relevant in fluid mechanics, non- linear optics, quantum gravity [4 , and SUSY QM 
applied to the harmonic oscillator [Sl[3]. In fact, a simple method to generate solutions of Pjy can be 
' supplied, since the SUSY technique provides explicit expressions for the relevant extremal states and 

^ . formulae connecting them with the Pjv solutions JTjiSj . 

This contribution is organized as follows: in Sections 2 and 3 we shall present the general framework 
of SUSY QM and polynomial Heisenberg algebras (PHA). In Section 4 we will generate the complex 
solutions to P/y, while in Section 5 we will study the eigenfunctions and eigenvalues of the non-hermitian 
Hamiltonians. Our conclusions shall be presented in Section 6. 

2 Higher-order SUSY QM 

In the fc-th order SUSY QM one starts from a given solvable Hamiltonian 



Ho^^l^ + Voi^), (1) 
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and generates a chain of standard (first-order) intertwining relations j9lfT0] 



Hence, the following equations must be satisfied 



l,...,fc. 



(2) 
(3) 



(4) 
(5) 



We are interested in the final Riccati solution afe(x, e^), which turns out to be determined either by fc 
solutions a\ (x, tj) of the initial Riccati equation or by fc solutions Uj of the Schrodinger equation obtained 
through the change a\(x, e^) — ul^juf 



3 General framework of PHA 



(6) 



A m-th order PHA is a deformation of the Heisenberg-Weyl algebra of kind [71[TT] : 

= ±L±, [L-,L+] = Qra+liH + 1) - Qm+liH) = PmiH), 

m+1 

Qrn + l{H) = L + L- = n " ' 



(7) 
(8) 



where Qm+iix) is a (m + l)-th order polynomial in a;, Pm{x) is a polynomial of order to in a; and Ei are 
the zeros of Qm+i{H), which correspond to the energies associated to the extremal states of H. 

For the second-order PHA (to = 2), let us suppose from now on that L+ is a third-order differential 
ladder operator of the form: 



(9) 



Using the standard first and second-order SUSY QM and decoupling the resulting system gives rise to 
equations for /, h, and V in terms of g and the following condition: 



2g 2 



2(x'^ - a) g+ -. 

9 



(10) 



This is the Painleve IV equation [Piv) with parameters a = £2 + ~ 2£i — 1 and b = — 2(£'2 — . 
Hence, if the three quantities £i are real, we will obtain real parameters a, h for the corresponding P/y. 

The key point now comes from the inverse method: instead of solving directly Piv to get a system 
ruled by a second-order PHA let us choose a system having this algebra (i.e. with third-order ladder 
operators) in order to obtain then the solutions to Pjv- 



4 Complex solutions to P/y with real parameters 

Let us use a theorem stating the conditions for the hermitian higher-order SUSY partners of the har- 
monic oscillator Hamiltonian to have second-order PHA (see [5j). The main requirement is that the k 
Schrodinger seed solutions have to be connected in the way 

Uj ^ [a^y^'^ui, ej=ei-(j-l), j^l,...,k, (11) 



2 



Complex solutions to Pjv through SUSY QM 



D. Bermudez, D. J. Fernandez C. 





20 






6 








15 










I (a 






s 


10 








— 2- /; 

/' \ 

H -^^"^^^^^^^ / ; 








is 


5 

























60 M 










-5 








-2- 1 










-10 






-4 









-4 -2 2 4 -4 -2 2 4 



X X 

Figure 1: Real (solid curve) and imaginary (dashed curve) parts of some complex solutions to P/y. The 
left plot corresponds to a = —5/2, b ~ —121/2 (fc = 2, ei = 5, A = 1, k = 5) and the right one to a = 9, 
63 = -2 (fc = 1, ei = 5, A = K = 1). 



where a~ is the standard annihilation operator, the only free seed ui is a nodeless real solution of Eq. ^ 
associated to a real factorization energy ei such that ei < Eq — 1/2. 

In order to overcome the last restriction, so that ei > Eq but the SUSY transformation being still 
non-singular, we need to use complex seed solutions. Let us choose them as complex linear combinations 
of the two standard linearly independent real solutions so that, up to an unessential factor, we will 
employ [T^ : 

u(x;e) = e-"'/' 

where A, k € M and iFi is the confluent hypergeometric (Kummer) function. 

Hence, through this formalism we will obtain simultaneously the fc-th order SUSY partner potential 
Vk{x) of the harmonic oscillator and the corresponding P/y solution g{x; ei), both of which are complex, 
in the following way 

Vkix) = xy2 - {ln[W^(ui, . . . , Uk)]}", g{x- ei) = -x - {H^e^ i^)]}' , (13) 

which can be compared with those obtained in [13] through Backlund transformations. Note that the 
extremal states of Hk and their corresponding energies are given by 

V-f, oc B+e-^'^\ #3 cx B+a+ui, (14) 

£2 = \, ^3 - £1 + 1, (15) 

where = . . . . The real and imaginary parts of the complex solutions g{x\ a, b) for two particular 
choices of real parameters a, b are plotted in Fig. [T] . 



iFi 



l-2e 1 



+ a;(A + iK) iFi 



3-2e 3 



(12) 



£1 - 



W{ui, . . . ,Uk-i) 
W(ui,...,uk) ' 

fife = £1 - (A: - 1), 



5 Non- hermit ian Hamiltonians 

Let us analyze the Hamiltonian Hk obtained by the complex SUSY transformation. The real case with 
ei < 1/2, leading to hermitian Hamiltonians, has been studied previously [9j[T4]. The spectrum Sp(iJfe) 
consists of an infinite ladder plus a finite one: there are two extremal states (both annihilated by L~) 
from which the two ladders start, one associated to and the other to Eq = 1/2; since the ladder starting 
from Ck ends at ei, the eigenfunction associated to ei is annihilated by L+. The actions of onto any 
other eigenstate of Hk are non-null and connect only eigenstates belonging to the same ladder. 
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Figure 2: Spectra of Hk and Hq for ei > 1/2 and (a) ei 7^ Ej; (b) ei — Ej. In (a), Sp{Hk) still contains 
one finite and one infinite ladders, but this is not true for (b). The dark bars represent the original and 
mapped eigenstates of Hq and Hk, while the light ones do for the k new levels of Hk- 

In the complex case, the new Hamiltonians necessarily have complex eigenfunctions (square-integrable), 
although the associated eigenvalues are still real. The resulting spectra for the non-hcrmitian Hamilto- 
nians Hk obey the same criteria as the real case, namely, they are composed of an infinite ladder plus a 
finite one which now could be placed, either fully or partially, above Eq — 1/2. 

The case when ei belongs to the spectrum of the original harmonic oscillator Hamiltonian is worth 
of a detailed study. Thus, let us consider a fc-order SUSY transformation with ei = Ej, j > k, and ui 
given by Eq. (|12p. i.e., ui is a complex linear combination of the eigenfunction -0^ of Ho and the other 
linearly independent solution of the Schrodinger equation. It is straightforward to see that the action of 
the third-order ladder operators applied to the eigenstate with energy Ei is given by 

L-{B+^i) cx B+0,_i, for / ^ {0,j - fc + 1}; L-{B+^o) = = 0, (16) 

L+{B+^Ji) cx B+in+i, for / ^ j; L+{B+^jj) = 0. (17) 

These results do not match with the established criteria for the non-singular real and complex cases with 
ei 7^ Ej since now it turns out that: 

E,+,, L-{B+^,+,)^B+il;,^0, (18) 
E,.k, i+(B,tV'j-fc) a B+^j_fc+i ^ 0, (19) 

where the shown energies correspond to the departure state. The resulting Hamiltonian is isospectral to 
the harmonic oscillator but with a special algebraic structure since now two states (those associated to 
Ej+i and Ej-k) are connected just in one way with the adjacent ones (see a diagram representing this 
in Fig. [2]). 

6 Conclusions 

We have described a simple method to obtain real and complex solutions g(x; a, b) of the Pjv with 
real parameters a, b. We have studied the properties of the resulting solutions, and we have analyzed 
the algebras, the eigenfunctions and the spectra of the non-hermitian SUSY generated Hamiltonians. 
F\irther extensions of this technique to obtain Pjv solutions associated to complex parameters a, b are 
under study. 
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